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Abstract 

(D : 

. We analyze the CP violating ratio e' je and rare K and B decays in the MSSM 

with minimal flavour and CP violation, including NLO QCD corrections and impos- 
^ . ing constraints on the supersymmetric parameters coming from e, i?^^ — B^^^ mix- 

. ings, B Xsj, Ag in the electroweak precision studies and from the lower bound on 

the neutral Higgs mass. We provide a compendium of phenomenologically relevant 
formulae in the MSSM. Denoting by T{Q) the MSSM prediction for a given quantity 
normalized to the Standard Model result we find the ranges: 0.53 < T{e' /e) < 1.07, 
0.65 < T{K+ TT+uu) < 1.02, 0.41 < r(i^L ^ 7r°i/z^) < 1.03, 0.48 < T{Ki^ 
7r°e+e") < 1.10, 0.73 < T{B Xsvv) < 1.34 and 0.68 < T{Bs fifi) < 1.53. 
We point out that these ranges will be considerably reduced when the lower bounds 
on the neutral Higgs mass and tan /? improve. Some contour plots illustrate the 
dependences of the quantities above on the relevant supersymmetric parameters. 
As a byproduct of this work we update our previous analysis of e' /e in the SM and 
find in NDR e'/e = (9.2t|g), a value 15% higher than in our 1999 analysis. 



1 Introduction 



One of the present central issues of particle physics phenomenology is the question whether 
the size of the experimentally observed direct CP violation in — > tttt decays can be 
described within the Standard Model (SM). 

The most recent experimental results for the ratio s'/s, which measures the size of the 
direct CP violation in Kl — > tttt relative to the indirect one, read: 



Together with the older NA31 measurement ((23 ± 7) • 10~^) [3], these data confidently 
estabhsh direct CP violation in nature and taking also the E731 result ((7.4 ± 5.9) • 10"^) 
[4] into account one finds the grand average 



close to the NA31 result. 

There are different opinions whether the grand average in (1.2) can be accommodated 
within the SM. We hst the results from various groups in table 1. The labels (MC) and 
(S) in the second column stand for two error estimates: "Monte Carlo" and "Scanning" 
respectively. The result from the Munich group given here is an update of the analysis 
in [5] done in the present paper. Similarly the result of the Rome group is the most 
recent estimate given in [6]. We do not list the values of the non-perturbative parameters 
i?6^^^'' and B^^"^^ used in different analyses as the comparison of these values could be 
misleading. The reason is that in certain papers these factors are independent and 
in other papers they depend on m^. Typically Bq'^'^^ = 0.8 — 1.6 and B^/^"^^ = 0.6 — 1.2. 
Exceptions are the analyses in [7] and [8] where Bq'^'^^ in the ballpark of 3 can be found. 

In [5, 6, 9] e'/e has been found to be typically a factor of 2 below the data and the 
KTeV result in (1.1) could only be accommodated if all relevant parameters were chosen 
simultaneously close to their extreme values. On the other hand the NA48 result is 
essentially compatible with [5, 6, 9] within experimental and theoretical errors. Higher 
values of e'/e than in [5, 6, 9], in the ballpark of (1.2), have been found in [7, 8, 10-13]. 
The result in [14] corresponds to Bq'^'^^ = B^^^'^^ = 1 (see section 7.1) as the Dubna-DESY 
group has no estimate of these non-perturbative parameters. Recent reviews can be found 
in [6, 10, 15]. Furthermore it has also been suggested that final state interactions (FBI) 
could enhance e'/e by a factor of two [10, 13, 16]. A critical analysis of this suggestion 
has been presented in [17]. In our opinion the issue of the size of FSI in the evaluation of 
e'/e is still an open question. An interesting recent proposal in [18], if realized, could put 
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(1.1) 



Re(£7£) = (19.2 ± 2.4) • 10"^ with xV^df =11.1/5 



(1.2) 
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Table 1: Results for e'/e in the SM in units of lO"''. 



the inclusion of FSI in the lattice calculations of hadronic matrix elements in principle 
under control. 

While all theoretical analyses quoted above use the NLO short distance Wilson co- 
efficients of [19-23], they differ in the evaluation of the hadronic matrix elements of the 
relevant operators. Being of non-perturbative origin, the latter calculations suffer from 
large theoretical uncertainties, which at present preclude any firm prediction of e'/e within 
the SM and its extentions. In this context, the matching between the short distance calcu- 
lations of Wilson coefficients and the long distance estimates of hadronic matrix elements 
is not fully under control. In particular, it is crucial that these two calculations are 
performed in the same renormalization scheme [24]. Moreover it has been pointed out 
recently in [25] that the presence of higher dimensional operators can in the case of low 
matching scales complicate further this issue. Finally, there is the question of isospin 
breaking effects that could be more important than initially thought but are difficult to 
estimate [26, 27]. There is a hope, however, that in the coming years the situation may 
improve through advanced lattice simulations and new analytical studies. 

In any case in view of the very large theoretical uncertainties and sizable experimental 
errors there is still a lot of room for non-standard contributions to e'/e. Indeed the results 
from NASI, KTeV and NA48 prompted several analyses of e'/e within various extensions 
of the SM like general supersymmetric models [28-32], models with anomalous gauge 
couplings [33] , four-generation models [34] and models with additional fermions and gauge 
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bosons [35]. Unfortunately several of these extensions have many free parameters and 
these studies are not very conclusive at present. The situation may change in the future 
when new data from high energy colliders will restrict the possible ranges of parameters 
involved. For instance a very recent analysis of the bounds on anomalous gauge couplings 
from LEP2 data indicates that substantial enhancements of e' jz from this sector are very 
unlikely [36]. 

On the other hand it appears that the e' jz data puts models in which there are new 
positive contributions to e and negative contributions to e' in difficulties. In particular 
as analyzed in [5] the two Higgs Doublet Model II can either be ruled out with improved 
hadronic matrix elements or a powerful lower bound on tan/3 can be obtained from e' je. 
In the Minimal Supersymmetric Standard Model (MSSM), in addition to charged Higgs 
exchanges in loop diagrams, also charginos contribute. For a suitable choice of the su- 
persymmetric parameters, the chargino contribution can enhance e' je with respect to the 
SM expectations [37]. Yet, as found in [37], the most conspicuous effect of minimal su- 
persymmetry is a depletion of e' je. The situation can be different in more general models 
in which there are more parameters than in the two Higgs doublet model II and in the 
MSSM, in particular new CP violating phases. As an example, in general supersymmetric 
models e' je can be considerably enhanced through the contributions of the chromomag- 
netic penguins [28-30, 38], Z^-penguins with the opposite sign to the one in the SM [30, 
39-41] and isospin breaking effects [31]. 

The purpose of the present paper is a new analysis within the MSSM with minimal 
flavour violation, which updates and generalizes the work of Gabrielli and Giudice [37]. 
Compared to the latter paper, the main new ingredients of our analysis are: 

• The inclusion of NLO-QCD and QED corrections; 

• the update in the CKM and non-perturbative parameters; 

• the imposition of constraints on supersymmetric parameters coming in particular 
from in the electroweak precision studies and from the lower bound on the 
neutral Higgs mass, in addition to the ones provided by the B decay rate; 

• the simultaneous study of the rare decays n^uu, Ki^ — >• 7r°z/z/, Ki^ vr^e+e", 
B Xsfi' and Bg /ip, which are sensitive to Z°-penguin contributions. 

As in [37] we will consider the MSSM with minimal flavour violation, in which the 
effective local operators are the same as in the SM and the only source of CP violation is 
the KM phase. Since the supersymmetric particles have masses C(Mw) or higher and the 
effective low energy operators are the same as in the SM, the impact of new contributions 
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in a renormalization group analysis is only felt in the initial conditions for the Wilson 
coefficients taken usually at /x = (9(M-w)- The renormalization group transformation 
from fi = (9(Mw) down to fi = 0{1 GeV) is on the other hand the same as in the SM. 
This simplifies the inclusion of QCD corrections considerably. 

In view of the fact that the new contributions enter only the initial conditions for the 
Wilson coefficients it is useful to follow [42] and to cast the expression for e'/e into a 
formula which exhibits transparently the new super symmetric contributions. One finds 
[5] 

e'/e - ImXt ■ F,, , F,, = Pq + ^ P,F, (1.3) 

r 

where = VtdV^*. The coefficients Pi depend on non-perturbative parameters and include 
NLO-QCD and QED corrections. They are common to the SM and the MSSM. The 
functions Fr resulting from various box and penguin diagrams contain both the SM and 
supersymmctric contributions involving new exchanges such as charged Higgs particles, 
charginos, squarks etc. They depend, in addition to nit, on the masses of the new particles 
as well as on a number of other supersymmctric parameters. The explicit expressions for 
Fr and the numerical values for the coefficients Pi will be given in sections 3 and 4 
respectively. 

There are many new contributions in the MSSM such as charged Higgs, chargino, neu- 
tralino and gluino contributions. However, in the case of minimal flavour and CP violation 
it is a good approximation to keep only charged Higgs and chargino contributions. 

Our paper is organized as follows. In section 2 we review the elements of the MSSM 
which are relevant for our analysis. In section 3 we give the list of all F^ functions which 
can be decomposed into W, charged-Higgs and chargino contributions 

Fr — {Fr)sM + {Fr)H + {Fr)x = {Fr)sM + {Fr)sUSY , (1-4) 

and we compare our results with those given in the literature. In section 4 we collect the 
formulae for e, B^^ — B^^ mixings, e'/e and the rare decays in question. In section 5 
some general considerations are given. In particular we address the issue of a universal 
unitarity triangle [43] for all models, like the SM and the constrained MSSM, that do 
not have any new phases beyond the CKM phase and no new operators. In section 6 we 
discuss constraints on super symmetric parameters coming from the B — > decay, 
in the electroweak precision studies and the lower bound on the neutral Higgs mass. The 
numerical analysis of the formulae of section 4 including the constraints from section 6 is 
presented in section 7. Finally section 8 contains a summary and conclusions. A list of 
auxiliary functions is given in the Appendix. 
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2 The Minimal Supersymmetric Standard Model 



2.1 FCNC Processes and CP Violation in the MSSM 

In general, the MSSM contains more than a hundred new parameters with respect to 
the SM. A large number of these parameters can be a new source of Flavour-Changing 
Neutral Currents (FCNC) and of CP violation. Phenomenological studies [38, 44-48] have 
shown that most of these new flavour- and CP-violating parameters arc constrained to be 
very small by the experimental results on FCNC and CP violation. One is then left with 
the so-called "flavour problem" and "CP problem" , i.e. with the difficulty of explaining 
the smallness of these new parameters. We do not discuss here the possible solutions 
to these problems, but rather assume that they have been solved leading to a version of 
the MSSM where no new flavour- and CP-violating parameters appear (minimal flavour 
violation). Although this is not a necessary outcome of the solution of the flavour and CP 
problems, it is certainly a well-defined and interesting scenario that leads to "minimal" 
modifications of the SM predictions on FCNC and CP violating processes. 



2.2 The Model 

Let us now specify the details of the MSSM with minimal fiavour and CP violation. 
After applying to the superfields those field redefinitions that diagonalize the quark mass 
matrices and M^, the squark mass matrices read 

^2 ^( Ml + {M^f + M| cos 2/5 ai") cot /? + At ) M„ \ 

" V + ^4 + (^")' + cos 2/5 a^n^ / 

^ / V^lV + {Maf + Ml cos 2/5 af (-//* tan (3 + A\) M, \ 

\ Md (-// tan /3 + A^) + {M^f + M| cos 2/3 ) 

where each entry denotes a three-by-three submatrix and where 

Au) _ 1 2 - 2 o Au) _ 2 - 2 p 

id) 1 , 1 • 2 m (d.) 1 • 2 m 

a}, =— 2 + 3sm Bw a)^ = — gSm Bw ■ 

The fact that the exphcit squark mass matrices in the two left-left sectors are related by 
the CKM matrix y is a remnant of SU{2) gauge symmetry. The 3x3 matrices M?^, 
M|^, Mj^, Au, Ad are a priori unconstrained complex matrices (except for the hermiticity 
of M|^, M?^, iu) P^^^^^t^^i^i^g soft supersymmetry breaking in the squark sector. 
We also define, as usual, tan P — V2/V1, while /i denotes the supersymmetric Higgs mixing 
parameter. 



(2.3) 
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In the basis chosen, all neutral-current couplings (couplings of neutral gauge bosons, 
neutral scalars, and neutral higgsinos and gauginos to quarks and squarks) are still diag- 
onal in the sense that they do not mix members of different quark supermultiplets, while 
the charged current couplings contain the CKM matrix in the same way as in the SM. 
However, FCNC arise after diagonalization of the squark mass matrices as soon as any 
one of the eight 3x3 submatrices of eqs. (2.1) and (2.2) is nondiagonal, while new CP 
violation effects are present if they are not real. Furthermore, we can only have and 
V^M~^V diagonal at the same time if 

M?^ ^rfiH. (2.4) 

Our assumption of minimal flavour and CP violation is consequently equivalent to real /x, 
real and diagonal A^, A^, M?^, Mj^, and = m^l. 

The elements of Au, A^ except the one multiplying the top quark mass are further 
constrained to be very small by requiring that the model does not develop charge- or 
color-breaking VEVs or potentials unbounded from below [49-51], so that we may neglect 
them for all squarks but the stops. Similarly, if tan /3 is not too large, the off-diagonal 
terms proportional to /i are negligible for q ^ i. In this case, the squark mass matrices 
are already diagonal except for mixing in the stop sector. The real, symmetric stop mass 
matrix is diagonalized via the orthogonal transformation 

~ o I m"? \ / cos 9f sin 9z\ , , 

TMfT^ = ' , , T - * * . (2.5 

* \ ml J' \-smei cosej ^ ' 

As a further simplifying assumption, we will set all other squark masses equal, which 
then implies rnl^ ^^ ~ for i — u,d, s, c, b. A strict equality is inconsistent with (2.1) 
and (2.2) for the "left-handed" squarks because of the different D-terms in the up and 
down sectors, but this difference is numerically insignificant and will be ignored here. We 
therefore set 

'^iL,iR = ^ = ^' ^' ^- (2-6) 

We then take the input parameters in the squark sector to be the two stop masses and 
the stop mixing angle %. Prom the above we also have 

= -rrif + m\ cos^ 9i -\- ml sin^ 9f (2.7) 

as a formula for computing m. We choose ^2 to be the light stop. 

In the leptonic sector we neglect Left- Right mixing, and assume, analogously to (2.4) 
and (2.6), 

m-^=m|, i^u,e. (2.8) 
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The chargino mass matrix is diagonalized by two orthogonal 2x2 matrices U and V, 
according to: 

U[ ^ ^ ^]v^=[ ^ , (2.9) 



rriwV^ cos P ji ) \ 



m 



where and tan (5 have been introduced before and M is the weak gaugino mass. M and 
jjL are taken to be real to avoid the appearance of new sources of CP violation. We prefer 
to choose the lightest chargino mass, which we call m^^, as an input parameter instead 
of M. 

In the Higgs sector, the only additional relevant parameter is the mass of the charged 
Higgs, mH+- 

Altogether, we are left with the seven input parameters Wf^, m^^, ^j, f^xii A*' tan/3, 
and mH+- The ranges for these parameters used in our numerical analysis arc discussed 
in section 7.2. In the computation of the rare decays n'^uu, Ki^ Tr^iyu, Ki^ 

Tr'^e'^e", B XgiyP and Bg ^ /ip, one more input is needed, the slepton mass m^. 

2.3 Chargino Interactions 

Let us denote by d the quarks {d, s, b) and by Uk the squarks {u, c, t) with k — 1,2. The 
interaction of charginos (x/, j = 1, 2) with d and Uk is then given explicitly as follows: 

= 9V:J {ZjkPL + Y^kPn) C {^^Uk + h.c. (2.10) 

where Pl,r = (1 =F 75)/2, while C is the charge conjugation matrix. The couplings Zjk 
and Yjk are real so that the only phases are present in the CKM elements. We have: 

v2M^smp 

777- ~ ~ 

^^■2 = /oa/ ■ R ^^^ ^" + ^^-1 ^" (2.12) 
V2Mvi/ sm (3 

= J"'^ / ,2Cosg, = (2.13) 
V 2Miy cos [3 

Z32 = - ^ ^.-2 sin = (2.14) 

V 2Mvi/ cos jj 

with V,j and f/jj defined in (2.9). The m„ denotes generally the mass of the up-quark 
(m, c, t) corresponding to a given up-squark u. 9^ are the mixing angles in the squark 
sector. Note that since we neglect quark masses other than that of the top, we have 
= for M 7^ t (see the preceding section) as well as Zj^ = 0. An exception is the 
B decay in which mi, and consequently also Zjk cannot be set to zero. 
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The Lagrangian (2.10) applies also to lepton-sneutrino-chargino interactions after the 
CKM factor has been removed and 6^, Yj2 and Zj^ all set to zero. If the masses of charged 
leptons are neglected only the couplings Yji = —Vji are non-vanishing in this case. 

Similarly the interaction of charginos with u and dk is given by 

= gVudU {Zj^Pl + Yj^Pr) X+4 + h.c. (2.15) 

where 

Yji = -Uji (2.16) 

Yj2^Zj,^Zj2^0. (2.17) 

The Lagrangian (2.15) describes, mutatis mutandis, the neutrino-slepton-chargino in- 
teractions once the CKM factor has been removed. 



3 Basic Functions 



We hst next the basic functions which enter the formulae for £, — mixing, e'/e and 
rare decays in section 4. The functions relevant for e and S° — B^ mixing are given as 
follows: 



SsM{c,t) = S{Xcw,Xtw) 



SH{c,t) ^ 



X. 



-L2{XcH, Xtui 1) + 



4tan4/3^^ ^"'^ ' tan^/? 
S^{t,t) = f{t,t) + f{c,c)-2f{t,c) 

f{c,t) — ^ Xt^.;^^yia^5^t^5^c,,^-4-t'3(3:^4x3) ^CftXj' ^XiXj) 



L'2iXcwi -^tW 1 -^Hw) LiiXcwi ^tWl t 



(3.18) 



<,3=1,2 
h,k=l,2 



where Xab = ?7T.^/m^. The functions S and -^1,2,3 are given in the appendix. These results 
agree with [44, 52, 53]. They agree also with [37] after the erratum has been taken into 
account. 

The functions which include the contributions from photon-penguins (£)), Z-penguins 
(C), boxes with external dd (-B^*^)), uu (-B^")), vD (S^'^^)), ee (S(^^)) and gluon-penguins 
{E), are given as follows: 



D = DsAiiXtv 



tan^/3 



(3.19) 



J=l,2 
fc=l,2 



C = CsM{xt^) + :^CH{x,,) + 2 [%%{5.,A*,cW(x,,^^.,x,,^^.) 



h,fc=l,2 



1 



+ - 



(3.20) 
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where 



A* 

^hk 



cos^ 9l k — h 

— sinOicosdi k^h 



(3.21) 



sin^ dx 



k^h^2 



i.j = l,2 
fe=l,2 



(3.22) 



fc=l,2 



(3.23) 



k=l,2 



5^''^ = 5sM(xtH^) + T E ^ji'^ii^wxAyji,yihB^x^(^kxr''-x,:^xix,) 



k=l,2 



(3.24) 



(3.25) 



E = £;5M(xt^) + :^EH{xt.) + 2 E fe^H^f.^xKtJ - (i^^)] ■ (3-26) 



3=1,2 
k=l,2 



The resuhs for the SM and charged Higgs contributions are the same as in [52] and 
agree with most papers in the hterature. The chargino contribution to the D function 
given here agrees with [44, 54] and also with the erratum in [37]. Our resuh for the 
chargino contribution to the C function agrees with [44, 54] and the erratum in [37]. The 
C functions given in these papers at first sight seem to differ from each other and from 
our resuh (3.20). However, after the summations over i,j,h,k have been performed and 
the c part subtracted, the resuhs of these papers reduce to our resuh. In particular the 
constant terms in [54] disappear. Note that the convention for 6^ in [54] differs from ours 
by a sign. 

Next our chargino contribution to B^'^^ differs from the result of [37] by an overall sign, 
while we find the chargino contribution to B^""^ to be a factor of 2 larger than given in 
the latter paper. Finally our result for the chargino contribution to S^^^^ agrees with [54] 
and differs by an overall sign from [44] . In view of the fact that these box contributions 
are small these differences are unimportant for the whole analysis. 
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4 Compendium of Phenomenological Formulae 



4.1 Basic Formula for e 

The indirect CP violation in X — > tttt is described by the well known parameter e. The 
general formula for e is given as follows 

where 

ReAo 

with Aq = ^ (7r7r)/=o) and AMk being the Kl-Ks mass difference. The off-diagonal 
element M12 in the neutral i^-meson mass matrix represents the K^-K^ mixing. It is given 
in the MSSM by 

M12 = ^^FlBKniKM^ [K%Fe, + Xl\2Fu + 2A:A,*r/3FcJ (4.3) 
where Aj = and 

Fee = Ssm{c, c) = XeW (4.4) 
Ftt = SsM{t, t) + Snit, t) + ^^(t, t) (4.5) 

Fet^SsM{c,t) + '^SHic,t) . (4.6) 

In what follows we neglect the charged Higgs contribution to Fee- Due to the assumed 
mass degeneracy for squarks in the first two generations, chargino exchanges contribute 
only to Ff^. In what follows we will use the NLO results of [55, 56] for the QCD factors 
Vi- 

r]i = 1.38 r]2 = 0.57 rjs = 0.47 . (4.7) 

as obtained in the SM. It should be stressed that in the leading logarithmic approximation 
(LO) r]2 is indeed the same for all three contributions in (4.5). At NLO however this is 
no longer true because then also gluon corrections to various box diagrams have to be 
included and these are different for the three contributions in (4.5). These corrections 
are unknown at present. Since in the SM the NLO corrections to ri2 mainly remove 
the renormalization scale ambiguities present in NLO, we think that in view of other 
uncertainties in the supersymmetry sector this procedure of using 772 at NLO is justified. 
In the case of F^ the QCD corrections to charged Higgs exchanges differ from the SM 
results already in the leading order and are given hy rj^ = 0.21 [52]. 
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Finally Bk is a well known non-perturbative parameter, Fx is the fT-meson decay 
constant and 171^ the ii'-meson mass. The last term in (4.1) constitutes at most a 2% 
correction to e and consequently can be neglected in view of other uncertainties, in par- 
ticular those connected with Bk- 

4.2 Basic Formula for — Mixing 

The B° — B^ mixing is usually described by (AM)^^^, the mass difference between the 
mass eigenstates in the B^ — B^ system and the B^ — B^ system, respectively. (AM)^^^ 
is given as follows: 

(AM),,, = ^VBmB,,^{BB,,FlJM^Fu\Vt^a,s)\' (4.8) 

where Fu is again the function given in (4.5). Bb is a non-perturbative parameter analo- 
gous to Bk, Fb is the B-decay constant and rjB the QCD factor given at NLO in the SM 
by 1]B = 0.55 [55, 57]. Supersymmetric corrections modify this factor first at the NLO 
level. As they are unknown at present we will use the SM value for rjB also in the MSSM. 

4.3 Basic Formula for e'/s 

Using the functions hsted in Section 3, the formula for e'/e of [5] generahzes to the 
supersymmetric case as follows: 

- = ImAt • F,, (4.9) 

e 

where 

F,,^Po + PxX + PyY + PzZ + PeE (4.10) 

with 

X = C-45("), Y = C-B^'^\ Z^C+^D. (4.11) 

The functions D, C, B^'^\ and E are listed in (3.19)-(3.26). Retaining only the SM 
contributions in the latter functions one obtains the formula (2.38) of [5]. 

The coefficients Pj are given in terms of the non-perturbative parameters Bq^^'^^ and 
Sg^^^"* and the strange quark mass ms(mc) as follows: 

p. =rf +rf i?6-hrf i?8 (4.12) 

where 



137 MeV 



Rq = B, 



(3/2) 



137 McV 



(4.13) 
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Ag^ = 290MeV 


MS 


= 340 McV 


MS 


= 390 MeV 


i 


0) 

' i 


6 

' i 


8 


0) 

' i 


' i 


(8) 
' i 


0) 

' i 


' i 


8 

' i 





-S.122 


10.905 


1.423 


-3.167 




12.409 


1.262 


-3.210 




14.152 


1.076 


A 


0.556 


0.019 





0.540 




0.023 





0.526 




0.027 





Y 


0.404 


0.080 





0.387 




0.088 





0.371 




0.097 





Z 


0.412 


-0.015 


-9.363 


0.474 




-0.017 


-10.186 


0.542 




-0.019 


-11.115 


E 


0.204 


-1.276 


0.409 


0.188 




-1.399 


0.459 


0.172 




-1.533 


0.515 





-3.097 


9.586 


1.045 


-3.141 


10.748 


0.867 


-3.183 


12.058 


0.666 



Table 2: Coefficients in the formula (4.12) for various in the NDR scheme. The last 

^ ' MS 

row gives the ro coefficients in the HV scheme. 

B\ and B\ parameterize the matrix elements of the dominant QCD-penguin (Qg) 
and the dominant electroweak penguin (Qg) operator respectively. The numerical values 
of rf\ rf^ and r'f^ for different values of A^ at u = nif. in the NDR renormalization 
scheme are given in table 2. This table differs from the one presented in [5] in that 
^n+n' = 0.25 has been replaced by Qri+rj' = 0.16 in accordance with the most recent 
calculation in [58] which gives ^Irj+r)' ~ 0.16 ± 0.03 (see however [26]). All other input 
parameters are as in [5] where further details connected with the coefficients in (4.12) and 
the parameters B^^"^^ and B^^"^^ can be found. 

4.4 Rare K and B Decays 

The well known expressions for rare K and B decays in the SM can easily be generalized 
to the corresponding expressions in the MSSM. We have [24] 



Br{K+ n+uP) = 1.54 • 10"^ • (imAt • X^-^^))^ + (RcAc ■ + ReXt ■ X 



(up) 



, (4.14) 



where Xc — (9.8 ± 1.4) • 10~^ stands for the internal charm contribution evaluated at NLO 
in the SM [59]. ReAc = — A(l — A^/2) with A = 0.22. The supersymmetric contributions 
reside only in X^''"^\ 



Br{Ki, ttVp) = 6.8 • 10"^ • (imAt • X^'^^))^ 



(4.15) 



Br{Kj^ ^ 7r°e+e-)dir = 6.3 • 10-^(ImAt)^(^7^^ + f,y) , 



(4.16) 
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where 



-4Z 



sin^ 



.y(ee) 



(4.17) 



sin i^iJ^i "^w 

with Pq = 3.05 ± 0.08 including NLO corrections [60] and sin^ ©w = 0.23. In (4.16) we 
have only given the directly CP violating contribution to this decay. 
Finally 



Br{B X,vV) = 1.5 • 10 



-5 1^*. 



(4.18) 



Br{Bs = 3.3 • 10" 



r{Bs) 
1.6ps 



n 2 



210 MeV 



\Vts\ 

0.040 



(4.19) 



where t{Bs) and Fb^ is the life-time and the decay constant of the Bg meson respectively. 
The functions X^""^^ and Y^^^'^ are given as follows 



— (J _ 4^b''^^\ y^^^^ — C — B^^^\ 



(4.20) 



with C, E^'^^") listed in (3.20), (3.24) and (3.25) respectively The function Z is 

defined in (4.11). 



5 General Considerations 

5.1 Strategy A: Universal Unitarity Triangle 

The SM, the MSSM and any model in which all flavour changing transitions are governed 
by the CKM matrix with no new phases beyond the CKM one and no new local oper- 
ators beyond those present in the SM, share a useful property [43]. Namely, the CKM 
parameters in these models extracted from a particular set of data are independent of the 
basic functions of section 3, they are universal in this class of models. Correspondingly 
there exists a universal unitarity triangle. In particular a determination of CKM param- 
eters and of this universal unitarity triangle without the knowledge of supersymmetric 
parameters is possible. 

Indeed using (4.8) one finds 



\Vtd\ _ ^ /rn^ 



(5.1) 



\Vts\ ^^mB,\{AM) 

That is the ratio k depends only on the measurable quantities (AM)^^^, thb^^ and the 
non-perturbative parameter ^. Now to an excellent accuracy [61]: 



Vtd\ = \Va\\Ru Rt = J{l-Qy + v^ (5.2) 
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\Vts\^\Vcb\{l-\x^ + Q\^) (5.3) 

where 

Q^Q{1-\X% (5.4) 

with A = 0.22, Q and n] being Wolfenstein parameters. We note next that through unitarity 
of the CKM matrix and the present lower bound on (AM)^ one has both in the SM and 
in the MSSM < ^ < 0.5. Consequently \Vts\ deviates from \Vcb\ by at most 2.5%. This 
means that to a very good accuracy the length of one side of the unitarity triangle is given 

by 

Rt^l (5.5) 

independently of supersymmetric parameters and mt- In order to complete the determina- 
tion of Q and f} one can use sin 2/3 extracted either from the CP asymmetry in Bd ipKs 
or from K nuP decays [62]. Both extractions are independent to an excellent accuracy 
of the supersymmetric parameters and mt. Once Rt and sin 2/3 have been determined in 
this manner, g and fj can be found through [43, 63] 

fj = ^^sin2/3-r_(sin2/3), g ^ 1 - fjr+{sm2p) (5.6) 

where r±{z) — (1 ± y/l — z^)/z. In general the calculation of g and fj from Rt and sin 2/3 
involves discrete ambiguities. As described in [43, 63] they can be resolved by using further 
information, for instance coming from jl^b/Kbl and e, so that eventually the solution (5.6) 
is singled out. 

As an alternative to sin 2/3 one could use the measurement of \/g^ + fj'^ by means of 
iKife/Kil but this strategy suffers from hadronic uncertainties in the extraction of |Kt6/Kb|- 

We observe that provided {AM)g has been measured and sin 2/3 extracted from CP 
asymmetry in Bd ipKs or K ^ nui? one can determine the "true" values of f] and g 
independently of supersymmetric parameters. Since A and \Vch\ = AX'^ are determined 
from tree level K and B decays they are insensitive to new physics as well. Thus the full 
CKM matrix can be determined in this manner. The corresponding universal unitarity 
triangle common to the SM and the restricted MSSM can be found directly from (5.6). 

Having the universal values of g and f] at hand one can calculate e, e'/e, (AM)^, (AM)^ 
and branching ratios for rare decays. As these quantities depend on supersymmetric 
parameters the results for the SM and the MSSM will generally differ from each other 
and one will be able to find out whether the SM or the MSSM, if any of these two models, 
is singled out by the data. In our opinion the use of the universal unitarity triangle is 
the most transparent strategy for distinction between models belonging to this class. It 
will certainly play an important role once (AM)^ has been measured and sin 2/3 extracted 
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either from the CP asymmetry in ipKs or from K -rvv decays. Other strategies 

for the determination of the universal unitarity triangle are discussed in [43] . 

5.2 Strategy B: Present 

As (AM)s and sin 2/3 from — > '4>Ks or X — > t^vV are presently not available we do 
not know the "true" values of the CKM parameters. However, by fitting theoretical 
expressions for e and (AM)^ to the experimental data we can find out for which values 
of CKM parameters the SM and the MSSM reproduce these data. As this determination 
involves the basic functions of section 3, the resulting values of CKM parameters in the 
MSSM will generally differ from the corresponding values in the SM. This difference has 
to be taken into account when making predictions for e' /e and rare decays branching 
ratios. 

A different strategy has been adopted in [64]. These authors have assumed all CKM 
parameters in the MSSM to be equal to those in the SM. This would be a correct strategy 
if we knew the universal triangle but in its absence it could be misleading. Indeed as a 
result of this assumption (AAf)^^^^ in [64] can differ from (AM)^^ by as much as 20% 
which is well outside the experimental error for {/\M)d that amounts to ±3%. The same 
comment applies to the parameter e where the experimental error is even smaller. This 
means that either the SM or the MSSM is already ruled out which of course is not true 
at present. 

In our opinion the only correct procedure at present is to impose the conditions 

{^M)f = (AM)^^s^ = (AM)^"P , = = . (5.7) 

As this conditions are used in our paper, our results for the rare decays branching ratios 
differ from the ones in [64]. 

Now, from (AM)^^ = (AM)^ssm foUows that 



(F,,)mssm ^^^^^ 



'^tdlMSSM — rtd\sM7Tr-s ZMSSM " 



where the last factor takes into account the small difference between QCD factors in the 
MSSM and the SM which appears at the two-loop level. We will set rf^'^^^^ = rj^^ in what 
follows. 

As (Ftt)MssM > (-Ftt)sM one has 

I ^td I MSSM < l^tdlsM ■ (5-9) 

Similarly, 

(ImAt)MSSM < (ImAt)sM, |(ReAt)MSSM| < |(ReAt)sM| (5.10) 
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when the constraints from (AM),^ and e are imposed. These hierarchies agree with the 
recent analysis in [65]. 

On the other hand \ Vts\ can only be modified by supersymmetric contributions through 
the change of q in (5.3). As this is an 0{\^) effect and < ^ < 0.5 from unitarity and 
\Vub/Vch\ we have within 2.5% 

l^tslMSSM = l^tsIsM = \Vcb\- (5-11) 
These results have the following consequences: 

• Quantities sensitive to \Vtd\-i ImA^ and ReA^ experience a suppression in the MSSM 
relative to the SM through the CKM parameters which may or may not be com- 
pensated by the loop effects relevant for a given process. In particular we have: 



)mssm 



Bt{Kl 



7rOe-e^ 



SM 



\dir 
mSSM 



(ImAt) 



2 

SM 



^MSSM 
^SM 



(ImAt)MSSM (^7A + l/7y)MSSM 

(ImAt)iM iy7A + y7v)sM 



(5.12) 
(5.13) 



ie'/e) 



MSSM 



(ImA 



t J MSSM 



^MSSM 



SM 



{£'/£)sM (ImAt)sM 
Moreover in the approximation of neglecting the charm contribution to 

one finds using (5.8) 

2 



Br{K+ 7r+z/z/)MSSM (Ft 



ttlSM 



"^MSSM 
^SM J 



Br{K+ 7r+z/z/)sM {Ftt)MSSM 
Of course, in our numerical analysis we keep the charm contribution. 



(5.14) 
(5.15) 



Quantities sensitive to \ Vts\ depend to a very good approximation on supersymmetric 
effects only through loop effects relevant for a given process. In particular we have 

Br{B XsVV)mssm 



Br{B XsUp) 



SM 



^MSSM 
VuP 
^SM 



(5.16) 



Br(B, 



MSSM 



BriBs 

• From (5.12) and (5.14) we have 



till) 



SM 



(^V - )aissm 



{e'/e) 



SM 



Br{K^ 



Vee 
-'MSSM 

-^SM J 



n 2 



.MSSM 



A'. 



SM 



\ Br{Ki^ ^ TT^uu) 



SM 



F!> 



SM 



■-MSSM 



(5.17) 



(5.18) 



These observations should be helpful in understanding the numerical results presented in 
section 7. It should also be remarked that the formulae (5.12)-(5.18) are also valid for 
strategy A after the ImAj and have been removed from these equations. 
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6 Constraints 



6.1 The B — ^ Xgj constraint 

The latest data [66] for the inclusive radiative decay rate of the B meson are known to 
constrain severely the accessible super symmetric parameter space (see for instance [67, 
68]). This decay mode is controlled by the magnetic penguin operator Qy, whose Wilson 
coefficient at the weak scale can be written as 

Cr^iMw) = + + AW, (6.19) 

where A^jlj, A^^\ and A^'' are the contributions from W, charged Higgs, and chargino 
loops, respectively. 

For a light charged Higgs mass, is large and always adds to the SM contribution, 
driving the branching ratio towards very high values. Therefore, only large chargino-stop 
contributions that interfere destructively with the charged Higgs ones can accommodate 
the present data [69]. On the other hand, unless the charged Higgs is very heavy, this 
requires the existence of at least one light chargino and one light stop and restricts con- 
siderably the allowed SUSY parameter space. In our analysis we therefore include this 
constraint using the latest theoretical and experimental information. 

A complete NLO analysis of 5 — > Xg^y is by now available in the SM [70, 71] and in 
the 2HDM [71, 72]. The references to earher LO calculations can be found in [24]. As a 
result of these efforts, the error on the theoretical prediction coming from uncalculated 
higher orders is estimated to be of the order of 5%, while large parametric uncertainties 
are induced by the semileptonic BR, used as a normalization, the bottom and charm 
masses, the CKM factor, and to a lesser extent by the top mass and as- In the case of 
the SM and of the type II 2HDM the overall theoretical error is generally around 10%. In 
the numerical analysis we fully implement all the NLO contributions to SM and 2HDM 
components including the latest refinements [73]. 

For what concerns the purely supersymmctric contributions, there exist no complete 
NLO calculation. However, the case considered in the present paper (minimal flavour 
violation) has been studied in [74, 75], where the Wilson coefficients have been expanded 
in inverse powers of the gluino and heavy squark masses. This approximation is designed 
to capture two sources of large NLO effects: i) if large cancellations between charged 
Higgs and chargino contribution occur at LO, they can be spoiled at NLO ; ii) the exis- 
tence of large splitting among supersymmetric particles can lead to large non- decoupling 
logarithms. 

In our numerical study we employ the results of [74] whenever they are applicable: we 
define a light mass scale Mi as the heavier between rht^ , fh^-^ , and a heavy mass scale Mh 
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as the minimum of nig, mt-^ and m; if the condition Mh > 2M/ is satisfied, the chargino 
contributions to the Wilson coefficients are evaluated at NLO, otherwise we use the LO 
expressions. Similarly, we employ the NLO expressions only if \6t\ < vr/lO. 

The estimate of the theoretical error is expected to compensate for this different treat- 
ment. We adopt a rather conservative approach: we scan over all the unphysical scales, 
using 2.4 < jib, Jib < 9.6 GeV and 40 < nw < 160 GeV (notation of [71]) and combine 
in quadrature the parametric uncertainties. The two extreme values of the BR obtained 
combining linearly these two uncertainties are then compared with the CLEO 95% C.L. 
BR^ > 2.0 X 10-^ and BR^ < 4.5 x 10"! 



6.2 Constraints from elect roweak precision tests 

The good agreement between the SM predictions and the clectroweak precision data 
favours those types of new physics for which contributions decouple from the precision 
observables. As we also consider light superpartners, we need to take into account the 
tight constraints on the supersymmetric spectrum emerging from this agreement. 

In the case of the MSSM, several recent and thorough analyses are available in the 
literature [67, 76]. The most relevant effect is due to the mass splitting of the superpart- 
ners, and in particular of the third generation squarks. Indeed, large splitting between 
TTibj^ and rhtj^ would induce large contributions to the electroweak p parameter of the same 
sign of the standard quadratic top quark term. This universal contribution enters the Z'^ 
boson couplings and the relation between Mw, and a and is therefore significantly 
constrained by present data. Following the assumptions of section 2, the contribution of 
the squarks to Ap = 1 — 1/p is given at one-loop order by the simple expression [77, 78] 

AV-^) = ^[-sin^^,cos^^,F(m?,m|) 

+ cos^ Oi F {ml , m? J + sin^ 9^ F (m| , mf J] , (6.20) 
with the function F given by 

F{x, y)^x + y log - . 

x-y y 

Notice that A^p^°) in (6.20) vanishes for degenerate squarks and when one of the squarks 
becomes very heavy it is proportional to the square of its mass. Of course, squark contri- 
butions are not the only supersymmetric corrections to Ap. However, they are the only 
potentially large ones, when the present direct exclusion bounds are taken into account, 
and the other contributions have in general the same sign [78]. 
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In the cases of the 1^-boson mass and of the effective weak mixing angle sin^ B^, for 
example, a doublet of heavy squarks would induce shifts proportional to its contribution 
to p, 

6Mw/Mw ^ 0.72 Ap; 6 sin^ 0^^ ^ -0.33 Ap. (6.21) 

As the relative experimental accuracy for these quantities is 5x10^"^ and 7 x 10^^, re- 
spectively, it is easy to realize that the present sensitivity to Ap is at the level of 1x10^^. 
Of course, a global fit provides a better way to explore this sensitivity. The model- 
independent analysis of [79], based on 1998 data, finds a value (3.7 ± 1.1) 10~^ for the 
parameter ei ~ Ap, where the central value is close to the SM result, but favours smaller 
vahies than in the SM. As Ap,^usy is positive, it is tightly constrained. The analysis of [80], 
on the other hand, includes also the information relative to the restricted Higgs sector 
of the MSSM and leads to Ap^usy = -0.000410 00^3 at the 2a level, which again shows 
no deviations from the SM. Another alternative formulation in terms of the parameter 
T aAp leads to very similar results (0 < Tsusy < 0.2 at 95% CL [67, 76]). We will 
therefore require 

AV°) < 1.5 X 10-^ 

which appears to be a conservative bound, also taking into account that QCD effects 
tend to enhance significantly the squark contributions to Ap [81], and that present data 
improve on the 1998 results employed in [79, 80], especially on Mw- 

6.3 Direct exclusion bound on 

The tree level MSSM relation between the mass of the lightest CP-even Higgs boson and 
Mz, Mh < Mz \ cos 2/? I < Mz-, is subject to very large radiative corrections (a summary of 
the most recent results with exhaustive references can be found in [82]), which relax it into 
a much looser bound < 135 GeV. The precise value of depends sensitively on the 
supersymmetric parameters and in particular on tan the mass of the pseudoscalar. Ma 
— which at lowest order is linked to the mass of the charged Higgs by = M\ + — 
and the third generation squarks, which give the most important loop contributions. 

At LEP-200 the lightest MSSM neutral Higgs can be produced mainly via the Higgs- 
strahlung process e+e^ ^ h Z and the pair production process e+e^ ^ h A. The exper- 
imental searches have so far set a lower bound on the mass of the light neutral Higgs of 
about 88 GeV [83]. As this already excludes large regions of the supersymmetric param- 
eter space, especially if tan j3 is low, it is important to include also this constraint in our 
analysis. To this end we employ approximate formulas based on [84] ; unfortunately, they 
cannot be used for Milr > 2Mq and do not cover part of our parameter space. There- 
fore, for MtLR > 2Mq we use the full one-loop expression to compute the Higgs mass 



19 




Figure 1: Contour plot in the rht^ —0t plane. The regions above the dash-dotted and dot- 
ted lines contain 95% and 100% of the generated points respectively, without taking the 
constraints discussed in section 6 into account. The regions above the dashed and contin- 
uous lines contain 95% and 100% of the allowed points respectively, once the constraints 
are taken into account. 

[85]. We notice that for large splitting between the stop masses new large logarithms 
enter the analysis which may induce relatively large QCD corrections and may need to 
be resummed. To the best of our knowledge, this case has not been considered in the 
literature. However we believe this does not affect our main conclusions. 

The relevance of the constraints discussed in this section on the SUSY parameter space 
can be seen, for example, in fig. 1, where we present the contour plot of the light stop 
mass versus the stop mixing angle in the ranges described in section 7.2, with and without 
the constraints discussed above. Clearly, a large portion of the region where one of the 
stops is very light is excluded by the present constraints. 
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7 Numerical Analysis 



7.1 Standard Model Parameters 

In order to make predictions for e' je and rare decays we need the values of ImAt and 
ReAj. They can be obtained from the standard analysis of the unitarity triangle which 
uses the data for iKttl, (AM)^ and (AM)^, where the last two measure the size 

oi B\^ — mixings. The relevant formulae for e and (AM)^^ are given in sections 4.1 
and 4.2 respectively. We also recall that the present lower bound on (AM)^ together with 
(AM)d puts the following constraint on the ratio |Vtd|/|Vt5|: 




(7.1) 



with the non-perturbative parameter ^ defined in (5.1). As discussed in section 5 the 
constraint in (7.1) is valid both in the SM and in the MSSM. 

The input parameters needed to perform the analysis of the unitarity triangle within 
the SM are: 

iKbl, |Kb|, "^c, rnt, Bk, \[WdFBa, (7.2) 

Their values are given in table 3, where rrit refers to the running current top quark mass 
defined aX. fi = mf"'"^. It corresponds to mf'''^ = 174.3 ± 5.1 GeV measured by CDF and 
DO. 

It should be emphasized that B^, \[^dFBi and ^ being obtained from the calculations 
of hadronic matrix elements in QCD are unaffected by the presence of super symmetry. 
The extraction of the remaining parameters in (7.2) from the data could in principle be 
affected by supersymmetric contributions but at least within MSSM such contributions 
are safely negligible. 

The additional input parameters required to perform the analysis of e' /e within the 
SM are 

^MS' "^slmcj, Bq , (7.3) 

Their values arc given in table 3, where the values for and rrisimc) correspond roughly 
to a,(Mz) = 0.119 ± 0.003, and m^{2 GeV) = (110 ± 20) McV respectively. 

(1/2) (3/2) 

Again the non-perturbative parameters Bq and B^ are unaffected by the presence 
of supersymmctry. The extraction of A^ and ms{mc) from the data could in principle be 
affected by supersymmetric contributions but at least within the MSSM such contributions 
are safely neghgible. 

In summary the values of the input parameters given in table 3 are common to the 
SM and MSSM. On the other hand the resulting values of ImAt, ReA^, e' /e and branching 
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Quantity 


Central 


Jjjrror 


\Vcb\ 


0.040 


±0.002 


\Vub\ 


3.56 ■ 10"^ 


±0.56 ■ 10"^ 


rric 


1.3 GcV 


0.1 GcV 


rrit 


165 GeV 


±5 GcV 


Br 


0.80 


±0.15 




200 MeV 


±40 MeV 




0.471 ps-i 


±0.016 ps-i 


(AM), 


> 14.6 ps-i 


95%C.L. 




l.lb 


±0.07 


MS 


340 MeV 


±50 MeV 


ms(mc) 


130 MeV 


±25 MeV 


r(V2) 

-°6 


1.0 


±0.3 




0.8 


±0.2 



Table 3: Collection of SM input parameters. 



ratios for rare decays will be generally different in these two models due to new super- 
symmetric contributions in the one-loop diagrams as summarized in sections 3 and 4. We 
follow here the strategy B of section 5 as the universal unitarity triangle is unknown at 
present. 

The values of the parameters in table 3 are precisely the ones used in [5] except for the 
lower bound on (AM), which has been improved since then [86]. The relevant references 
to table 3 can be found in [5] . 

The only controversial values in table 3 are the ones for B^^'^\ B^^"^^ and mg. Other 
authors would possibly use different values for these parameters. As one of the main 
purposes of the present paper is to investigate the impact of supersymmetric contributions 
on our e'/e analysis within the SM [5], we will keep the values of all the SM input 
parameters as in [5] . Reviews of the values of B^^^"^^ , B^^"^^ and rris used by other authors 
can be found in [5, 15]. Wc will briefly comment at the end of this section how our analysis 
would change had we used different values for Bq^^"^^ and B^^'^\ 

7.2 Supersymmetric Parameters 

As discussed in section 2.2, we will choose the following supersymmetric input parameters: 
mti, mt2, % m^i, tan/3, m^ji, ji, mi, (7.4) 
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where mt^ and rfit^ are the heavy and Ught stop mass respectively, 61 is the mixing angle 
in the stop mass matrix, the light chargino mass, tan/? the ratio of the two Higgs 
vacuum expectation values, the charged Higgs mass, the Higgs superfield mixing 

parameter and is the common slepton mass entering the expressions for rare K and B 
decays. We now discuss the ranges in which we randomly generate with flat distributions 
the parameters in eq. (7.4). 

The region of SUSY parameter space in which we expect SUSY contributions to e'/ e 
to be maximal (in absolute value) corresponds to the case of a light stop and chargino. 
We therefore allow for the lighter stop and chargino to be as hght as allowed by present 
constraints, 



Concerning the stop mixing angle, we generate it in the range — 7r/2 < 6]; < 7r/2. The 



charge and colour breaking minima in the potential. To avoid this, following ref. [87] we 
impose the bound 



which for a typical Higgs mass Mh ~ 90 GeV amounts to m? > - (90 GeV)^. We also 
impose rh > 250: Since minimal flavour violation requires the soft breaking left-handed 
squark mass to be universal, this imphes rht^ ^ 300 GeV for 9^ ~ ±7r/2. 

We vary the charged Higgs mass in the range 100 GeV < m//± < 1 TeV, and choose 

< 500 GeV. In this work, we do not analyze possible effects in the large tan/5 regime, 
where additional operators can arise. We can therefore limit ourselves to the region 
1.2 < tan /3 < 6, since the dependence on tan P becomes very weak already for tan /3 > 3, 
until one enters the large tan P regime. Finally, we impose < 1 TeV and 90 GeV < 
rfii < 500 GeV. We do not enforce GUT relations between SUSY parameters since we are 
considering a more general class of theories described in terms of the SUSY parameters 
at the electroweak scale, with no reference to the structure of the theory at higher scales. 
The NLO calculation adopted for BR(b 57) introduces a dependence also on the gluino 
mass rhg, which we choose to vary in the range 250 GeV < rhg < 1 TeV. 

We scan the parameter space specified above and impose the constraints discussed in 
section 6. Then for the points allowed by these constraints we perform the determination 
of ImAt and ReAt from AF = 2 processes and then compute e'/e and branching ratios for 
rare decays. We provide the details of this computation below. 

The ranges for supersymmetric input parameters are summarized in table 4. 



90 GeV < rh^^ < 250 GeV, 90 GeV < rht^ < 800 GeV. 



(7.5) 



case of a very light and right-handed t2 requires m?^ < 0. This might however induce 




(7.6) 
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Quantity 


TV ,T • 

Mm 


TV /T 

Max 


rhti 


90 GeV 


1 TeV 




90 GeV 


800 GeV 


Oi 


-n/2 


7r/2 


rnxi 


90 GeV 


250 GcV 


rfig 


250 GeV 


i lev 


tan/3 


1.2 


6 


mH± 


100 GeV 


ITeV 


1^ 


-500 GeV 


500 GeV 


me 


90 GeV 


500 GeV 



Table 4: Collection of supersymmetric input parameters. 
7.3 Strategy 

The uncertainties in the input parameters of table 3 screen considerably the effects of 
supersymmetric contributions. In order to study various trends and patterns of super- 
symmetric effects in ImAj, ReAj, e' /e and rare decays we proceed as follows: 
Step 1: 

We set the SM parameters to their central values in table 3 and, using a random 
generation of ~ 6 • 10® points, calculate the ranges for the ratios 

^^^^^ ^ (^r)MSSM ^ ^ ^ ^{yV)^Y^ee) ^77) 

yFr)sM 

with (Fr)sM and (Fr-)MSSM representing the SM and total contributions to the functions 
in question respectively. S = with Fft given in (4.5). These ratios are of interest for 
both strategies of section 5. In order to investigate the impact of various constraints on 
SUSY parameters we consider the following cases: 

1. Only the constraints from e, B^^ — B^^- mixings and direct bounds on masses of 
squarks, sleptons, charginos and charged Higgs particles are imposed. 

2. The additional constraint from B — > Xg^ is imposed. 

3. The additional constraint from Aq is imposed without any constraint from B — > 

4. The additional constraint from the lower bound on the neutral Higgs mass is imposed 
without any constraint from B — > Xg^ and A^. 
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5. All constraints are taken into account. 

This step allows us to identify the dominant supersymmetric contributions. 
Step 2: 

We repeat the analysis of Step 1 calculating this time the ranges for the ratios 

rp,!!-. (g7g)MSSM 

^(^AJ^^^TT^n (7-8) 

and analogous ratios for ImA^, Re At, and the branching ratios for rare decays. 

We also investigate for which sets of supersymmetric parameters the maximal and 
minimal values of these ratios are obtained. 

Step 3: 

We repeat the analysis of Step 2 in the case of e'/e, investigating the dependence of 
T[e'lE) on B'il''\ 

7.4 Results of Step 1 

Our analysis shows that, taking only the constraints (1) into account, supersymmetric 

contributions to the box functions B^"^^ and B^^^^ amount respectively to at most 

4%, 7% and 11% of their SM values and are suppressed below 3%, 6% and 11% after 
all constraints are taken into account. These effects are only relevant in the case of the 
function X where the box function enters with an additional factor 4. The supersymmetric 
contributions to the functions li?*^^^)!, C, E and S can be as high as +34%, +51%, 
+93%, +107% and +159% respectively if no constraints from B ^^7, I\q and (Mh)min 
are imposed. While S is always enhanced, C, l^l and E can be suppressed up to 24%, 
7% and 17% respectively. The constraints (2-5) turn out to have only a minor impact on 
the functions B^'^^\ D, E and S. On the other hand the maximal enhancement of C is 
reduced to 29%. 

The function E is absent in n^uu, Ki^ n^uu, B Ksfv and Bs fJ^fi, and 

plays only a minor role in e'/e and a negligible role in Ki^ TT'^e~^e~. Consequently its 
modification is immaterial for our analysis. On the other hand the enhancement of S can 
have considerable impact on the values of ImAj and ReAj as we will see below. Similarly 
the supersymmetric effects in the functions C and D are sizable, although these effects 
are less visible in the functions X, Y, Z due to various cancellations between box and 
penguin contributions. 

In table 5 we show the minimal and maximal values of the ratios K{Fr) defined in (7.7) 
in the case of no constraints from B Xg^, Ag and {Mh)min and after the imposition of 
these constraints. The anatomy of these three constraints is shown in table 6. In the case 
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No Constraints 


All Constraints 


K{Fr) 


min 


max 


min 


max 




0.96 


1.00 


0.97 


1.00 




0.93 


1.01 


0.94 


1.01 




0.89 


1.01 


0.89 


1.01 




0.94 


1.34 


0.96 


1.34 


K{C) 


0.76 


1.51 


0.76 


1.29 


K[D) 


0.93 


1.93 


0.95 


1.87 


K{X) 


0.87 


1.26 


0.87 


1.15 


K(Y) 


0.81 


1.41 


0.81 


1.23 


T 7- f r-7\ 

K{Z) 


0.68 


1.65 


0.70 


1.35 


K{E) 


0.83 


2.07 


0.83 


1.85 


K{S) 


1.01 


2.59 


1.01 


2.13 




0.83 


1.25 


0.85 


1.15 




0.82 


1.41 


0.82 


1.24 




0.59 


1.27 


0.76 


1.23 



Table 5: The minimal and maximal values of the the ratios K{Fr) without constraints 
and with all constraints taken into account. 

of the functions E and 5" we observe the features just discussed. The supersymmetric 
effects in the functions X, Y and Z can be at most ±15%, ±23% and ±35% respectively 
after all constraints have been imposed. This is related to the fact that while the function 
C is positive, the function D is negative. Consequently the observed enhancement of |D| 
results in a suppression of Z. This suppression is reduced by the factor 1/4 with which 
the function D enters the function Z and by the enhancement of C, which for certain 
values of supersymmetric parameters can overcompensate the suppression of Z through 
D so that an enhancement of Z is possible. 

It should be stressed that in the two Higgs doublet model the function Z is always 
enhanced through charged Higgs effects. The possible suppression of Z in the MSSM 
comes then from chargino exchanges that are particularly important in the photon penguin 
function D. 

Let us now identify the region of SUSY parameter space where the maximal deviations 
from the SM arc possible. The function S takes its largest values for a light, mainly right- 
handed stop, with a large splitting between the stops, and small tan/3 (see fig. 2). The 
function C reaches its maximal values for small tan (3 and large stop mixing, while the 
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B — > 


only 


Af) only 


{Mh)nnn Only 


K{Fr) 


min 


max 


min 


max 


min 


max 




0.96 


1.00 


0.96 


1.00 


0.96 


1.00 




0.93 


1.01 


0.94 


1.01 


0.93 


1.01 




0.89 


1.01 


0.89 


1.01 


0.89 


1.01 




0.95 


1.34 


0.94 


1.34 


0.95 


1.34 


K{C) 


0.76 


1.51 


0.76 


1.39 


0.76 


1.51 


K{D) 


0.93 


1.93 


0.93 


1.93 


0.94 


1.87 


K{X) 


0.87 


1.26 


0.87 


1.20 


0.87 


1.26 


K{Y) 


0.81 


1.41 


0.81 


1.31 


0.81 


1.41 


/ r7\ 

K{Z) 


0.68 


1.65 


0.68 


1.47 


0.70 


1.65 


K[E) 


0.83 


Z.07 


0.83 


2.07 


0.83 


1.85 


K{S) 


1.01 


2.59 


1.01 


2.59 


1.01 


2.18 




0.83 


1.25 


0.83 


1.20 


0.85 


1.25 




0.82 


1.41 


0.82 


1.31 


0.82 


1.41 




0.61 


1.27 


0.68 


1.27 


0.59 


1.23 



Table 6: The anatomy of various constraints on the minimal and maximal values of the 
the ratios K{Fr). 

minimal values are attained for a hght, mainly right-handed stop and light chargino (see 
fig. 3). The function |D| is almost always larger than in the SM, and reaches its maximum 
for a hght, mainly right-handed stop, a hght chargino and small (see fig. 4). Therefore, Z 
is enhanced with respect to the SM for small tan (3 and large stop mixing, and suppressed 
for a hght, mainly right-handed stop and hght chargino (See fig. 5). 

As seen in table 6 the constraints from Ag and the lower bound on mass are 
slightly more effective than the constraint from B — > ^^7. Yet the comparison of tables 5 
and 6 shows that the inclusion of all bounds simultaneously eliminates a range of possible 
values for X, Y and Z which are allowed if each constraint is imposed separately. 

To evaluate the effect of the bounds discussed in section 6, we show in fig. 6 the ranges 
for K{S) and K{C) without applying the constraints: the comparison with figs. 2 and 3 
clearly shows the impact of the bounds. 

From the above considerations it is clear that any future increase in the lower bounds 
on Mh and tan/? will considerably reduce the possible enhancements of S and Z, and 
therefore, as we shall see in the following, the possible suppression of e'/e. 
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Figure 2: Contour plots of K{S) as a function of rfit^, rfit^, tan/3 and Of, including all 
the constraints discussed in section 6. The dashed and continuous lines contain 95% and 
100% of the allowed points respectively. 



7.5 Results of Step 2 

In table 7 we show the minimal and maximal values of the ratios T defined in (7.8) in 
the case of no constraints from B — > Xg^, Ag and {Mfi)min ^-nd after the imposition of 
these bounds. The anatomy of these three constraints is shown in table 8. We make the 
following observations, taking all the constraints into account: 

• The extracted ImA^ and |ReA(| are suppressed by supersymmetric contributions. 
This is related to the enhancement of Fu in the MSSM. 

• e'/e can be enhanced by at most 7% but can be suppressed even by a factor of 
2. This pattern agrees with the one found in [37] although the supersymmetric 
effects in e'/e found here are smaller than in [37] due to the tighter constraints on 
supersymmetric parameters imposed in the present analysis. We will elaborate on 
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Figure 3: Contour plots of i^(C) as a function of mtj, mtj, tan/3 and ^f, including all 
the constraints discussed in section 6. The dashed and continuous lines contain 95% and 
100% of the allowed points respectively. 



these results below. 

• The branching ratios for — > tt+i/P, — 7r°i/P and — > 'K^e^e" are basically 
suppressed by supersymmetric contributions relative to the SM expectations. The 
suppression of Br{K^ — > 'k'^vv) can be at most 35%. The corresponding suppres- 
sions in — > 'K^vv and Ki^ — > vr'^e+e" can be as high as a factor of 2.5 and 2 
respectively. As seen in (5.12), (5.13) and (5.15) and tables 5 and 7, these suppres- 
sions are related to the suppressions of ImAt and jReA^I that cannot be compensated 
by possible enhancements of the functions X, Y and Z. 

• As expected on the basis of (5.16) and (5.17) the pattern of supersymmetric effects is 
different for the rare decays B — > XgUP and Bg nfi. Here possible enhancements 
and suppressions can be directly calculated using (5.16), (5.17) and table 5. The 
supersymmetric effects in Br{B XgUu) can be at most ±34%. Br{Bs — > /i/2) can 
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Figure 4: Contour plots of K{D) as a function of fht^i "rh^^, /i, rht^, tan P and 9f, including 
all the constraints discussed in section 6. The dashed and continuous hues contain 95% 
and 100% of the allowed points respectively. 

be enhanced up to 53% and suppressed up to 32%. 
7.6 General Patterns: e'/e 

Let us elaborate on our results for e'/e. The dominant supersymmetric contributions to 
e and {AM)^^^ originate in the functions Snitjt) and S^{t,t) contributing to Fu in (4.5). 
These two functions add positive contributions to Fu which is positive in the SM. As a 
result, the values of ImXt and |ReAt| in the MSSM are smaller than in the SM for any 
choice of the supersymmetric parameters in (7.4). This is clearly seen in table 7. The 
suppressions of ImAj and |ReAt| are largest for a very light, mainly right-handed stop, for 
a large splitting between the stops, for small m^^ and small tan/3. 

The dominant supersymmetric contributions to F^/ in (4.10) originate in supersym- 
metric contributions to the function Z in (4.11) which results from Z'^-penguin diagrams 
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Figure 5: Contour plots of K{Z) as a function of rht^, rhti, tan/3 and 9f, including all 
the constraints discussed in section 6. The dashed and continuous lines contain 95% and 
100% of the allowed points respectively. 



and photon-penguin diagrams. The corresponding contributions to the functions X , Y, 
and E have a considerably smaller impact on F^/. This is related to a large extent to the 
fact that the coefficient Pz in (4.10) is substantially larger than the coefficients Px, Py 
and Pe- 

Setting the parameters of tables 2 and 3 to their central values we find 

F^/ = 10.3 -7.7- Z + 1.1 (7.9) 

where the last term represents the contributions of X , Y and E evaluated within the 
SM. The first term represents the part unaffected by supersymmetric contributions. It is 
dominated by the contribution of the QCD-penguin operators, whose Wilson coefficients 
are strongly enhanced by renormalization group evolution from scales C(Mw) to scales 
0{1 GeV). As stated above, the supersymmetric contributions affect F^i dominantly 
through the function Z, which equals 0.66 in the SM. 
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Figure 6: Contour plots of K{S) and K{C) as a function of iht^ and ^f, without imposing 
the constraints discussed in section 6. The dashed and continuous hnes contain 95% and 
100% of the generated points respectively. 

The charged Higgs contribution to the function Z is positive for any value of tan/3 
and mH±- Consequently, as seen in (7.9), this contribution suppresses F^i and e' /e. 
This suppression of e' /e is enhanced by the suppression of ImAj through charged Higgs 
contributions as discussed above. A detailed numerical analysis of these effects within the 
two Higgs doublet model has been presented in [5, 52] 

As pointed out in [37], depending on the choice of the supersymmetric parameters, the 
chargino contribution to the function Z and more generally to F^i can have either sign. 
Indeed as seen in table 5 the function Z can be suppressed up to 30% relative to its SM 
value even if all constraints are imposed. If this happens the suppression of e' /e due to 
Z°-penguins observed in the SM and in the two-Higgs doublet model becomes smaller. 

Whether e' /e can be enhanced in MSSM over the SM prediction depends on whether 
the positive contribution of charginos to the function F^i can overcompensate the corre- 
sponding negative charged Higgs contribution to this function and the suppression of ImA^ 
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No Constraints 


All Constraints 


T 


min 


max 


min 


max 


T(ImAt) 


0.57 


1.00 


0.66 


1.00 


T(ReAt) 


0.78 


1.00 


0.81 


1.00 


T{e'/e) 


0.42 


1.07 


0.53 


1.07 




0.59 


1.09 


0.65 


1.02 


T{Ki^ ^ ttVz/) 


0.28 


1.12 


0.41 


1.03 


T(Xl 7r°e+e-) 


0.33 


1.10 


0.48 


1.10 




0.70 


1.60 


0.73 


1.34 


T(5s ^ /x/i) 


0.68 


1.99 


0.68 


1.53 



Table 7: The minimal and maximal values of the the ratios T without constraints and 
with all constraints taken into account. 





B 


only 


only 




nm only 


T 


min 


max 


min 


max 


min 


max 


T{lmXt) 


0.57 


1.00 


0.57 


1.00 


0.64 


1.00 


T{ReXt) 


0.78 


1.00 


0.78 


1.00 


0.80 


1.00 


ne'/e) 


0.42 


1.07 


0.46 


1.07 


0.42 


1.07 




0.59 


1.08 


0.59 


1.02 


0.65 


1.09 




0.28 


1.12 


0.28 


1.03 


0.41 


1.12 


T{Ki^ 7r°e+e-) 


0.33 


1.10 


0.33 


1.10 


0.48 


1.10 


T{B Xsuiy) 


0.70 


1.60 


0.70 


1.46 


0.73 


1.60 


T{Bs /i/2) 


0.68 


1.99 


0.68 


1.73 


0.68 


1.99 



Table 8: The anatomy of various constraints on the minimal and maximal values of the 
the ratios T. 

due to supersymmetric contributions. While the chargino contribution to F^i can indeed 
overcompensate the corresponding negative charged Higgs contribution (see Table 5) the 
resulting enhancement of F^i is generally insufficient to cancel the suppression of ImAj. 
Consequently for the dominant part of the allowed values of supersymmetric parameters 
e' /e in MSSM is suppressed with respect to its SM estimate. This is clearly seen in the 
contour plots in fig. 7 where we show the range for e' / £ as a function of iht^ , rht^ , tan /3 
and with other parameters varied in the allowed ranges. The strongest suppression of 
e'/e takes place for lowest tan/3, mH±, for a light, mainly right-handed stop and for large 
splitting between the stops. The enhancement of e'/e over its SM value is found for large 



33 




1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 -0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 



tan(|3) e,/7t 

Figure 7: Contour plots of s'/s as a function of m^j, tan/3 and 6i, imposing the 
constraints discussed in section 6. The dashed and continuous curves contain 95% and 
100% of the allowed points respectively. The dashed horizontal hne corresponds to the 
SM prediction for central values of the parameters in table 3. 



tan/5 and and small m^^. In this case the charged Higgs contribution becomes much 
smaller than the chargino contribution and the enhancement of F^i through the chargino 
contribution can overcompensate the suppression of ImA^. 

Finally, let us comment on the dependence of the above results on the hadronic pa- 
rameters entering e' /e. It is well known that for central values of the parameters in table 
3 there is a strong cancellation between QCD and electroweak penguin contributions, as 
shown explicitly in eq. (7.9). In this case, the enhancement or suppression of electroweak 
penguins due to SUSY contributions has a relatively large impact on the value of e' /e. If, 
however, b'^^'^^ is larger than one, e' /e becomes more QCD-penguin dominated, the can- 
cellation becomes much less effective and therefore also SUSY contributions have a smaller 
impact on the value of e' /e (as discussed above, the main contribution to QCD penguins 
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No Constraints 


All Constraints 


d(1/2) 
^6 


T, K 


min 


max 


min 


max 


1.0 




0.59 


1.27 


0.76 


1.23 


i.U 


T[e /e) 


0.42 


1.07 


0.53 


1.07 


1.5 




0.78 


1.13 


0.87 


1.11 


1.5 


T{e'/e) 


0.53 


1.02 


0.60 


1.02 


2.0 


K{F,,) 


0.85 


1.08 


0.91 


1.07 


2.0 


Tie'/e) 


0.55 


1.00 


0.63 


1.00 



Table 9: The minimal and maximal values of the the ratios K{F^i) and T{e'/e) for different 
values of Bq'^'^\ without constraints and with all constraints taken into account. 



at the hadronic scale comes from the mixing with current-current operators in the RGE 
evolution, and is therefore independent of SUSY contributions). This is clearly displayed 
in table 9, where we report the extremal values of K{F^/) and T{e'/e) for Bq^^"^^ — 1, 
1.5 and 2. Similarly, increasing and decreasing the value of B^^^\ the effects of SUSY 
contributions to e'/e are increased and decreased respectively. 

7.7 Global Analysis 

In principle we could next vary all the SM parameters of tables 2 and 3 in order to obtain 
the full ranges for e'/e and for the branching ratios of rare decays in the MSSM. We do 
not think such an analysis would be very instructive at present. Roughly speaking the 
full ranges for the quantities considered in this paper can be obtained by using the SM 
ranges, given for instance in [24], and rescaling them by the enhancement and suppression 
factors of table 7. 

On the other hand we have updated our 1999 analysis of {e'/e)sM [5] by modifying 
simply = 0.25 to = 0.16 recently obtained in [58]. The result in NDR is given 

in table 1, to be compared with e'/e = of our previous analysis. 

8 Summary 

We have analyzed the CP violating ratio e'/e and rare K and B decays in the MSSM with 
minimal flavour and CP violation. A compendium of phenomenological formulae that 
should be useful for future analyses is given in section 4. In this version of the MSSM 
no new operators beyond those present in the SM contribute and all flavour changing 
transitions are governed by the CKM matrix. In particular there are no new phases beyond 
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the CKM phase. This imphes that the supersymmetric effects enter various quantities of 
interest exclusively through the one-loop functions. 

In spite of this simplification the analysis is complicated by the fact that in addition to 
a number of poorly constrained supersymmetric parameters one has to deal with hadronic 
uncertainties. The latter plague in particular e' je and to a lesser extent e and g — B^ ^ 
mixings. On top of this there are uncertainties in the CKM parameters. Therefore in order 
to uncover supersymmetric effects we have set the input parameters like \Vub\, \Vcb\, B^, 
b'^^'^\ Bf^'^\ that are insensitive to supersymmetric parameters, to their central values. 
Imposing subsequently constraints on the supersymmetric parameters coming from e, 
B^ g — B^ g mixings, B — > Xg^, Ag in the electroweak precision studies and from the lower 
bound on the neutral Higgs mass, we have calculated the MSSM predictions for e'/e and 
rare branching ratios normalized to the SM results. Our findings are summarized in table 
7 and in the abstract. The calculated quantities in the K-system are generally suppressed 
by supersymmetric effects, although a 7% enhancement of e'/e is still possible. In the 
case of S — > XgiyP and Bg ^ /ip, both suppressions and enhancements in the ballpark of 
30 — 50% are still possible. We have pointed out that an improved lower bound on the 
neutral Higgs mass could considerably reduce the allowed ranges. 

In view of very large hadronic uncertainties in e'/ e it will be very difficult to distinguish 
the MSSM prediction for this ratio from the one in the SM. If future improved calculations 
of Bq'^'^^ and B^g^'^^ will demonstrate indisputably that {e' /e)sM is below the experimental 
data our results show that the MSSM will be of little help to remedy this difficulty. On 
the other hand if it turns out that {e'/e)sM is too large then supersymmetric effects in 
the MSSM could help to obtain the agreement with data. 

As rare decays such as n^uu, K]^ n^uu, B — * X^z/z/ and Bg ^Ji are 

theoretically very clean, they offer a much better place to look for supersymmetric effects 
within the MSSM than e'/e. The analyses of these decays will be particularly transparent 
as soon as {AM)d/ (AM)s and the CP violating asymmetry in B ^ ipKs will be precisely 
measured allowing the construction of the universal unitarity triangle [43]. If it turns 
out that the data for ti^vv and Ki^ ti^vv are above the SM expectations, it 

will become clear than more complicated supersymmetric extentions of the SM than the 
constrained MSSM considered here are required to understand the experimental results. 
In the case of 5 — > XgVD and Bg /ip, no such clear conclusions will be possible as both 
enhancements and suppressions of these branching ratios are allowed. On the other hand 
with precise data the MSSM predictions could in principle be distinguished from the SM 
ones. 
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Appendix 



Here we collect the functions which enter the basic functions of section 3. We use the 
notation of [37, 52]. 

3 1 



S{x,y 

Li{x,y,z 
L2{x,y,z 
L3{x,y,z 

F{x,y,z 



xy 



1 3 1 

A^2l-x ~A{l-xY 



3 1 1 



\ogx 

X 



xy [F{x, y, z) + F{y, z, x) + F{z, x, y)] 
xy [xF{x, y, z) + yF{y, z, x) + zF{z, x, y)] 

—L2{x,y,z) 
xy 

xlogx 



{x — l){x — y){x — z) 



Box(A S=l) 



BUx,y,z) = — -Li{x,y,z) 
^ xy 

Bi{x,y,z) = L3{x,y,z) 



7-penguin 



^ , , a;2 (25 - 19a;) (-3x^ + 30x^ - Mx'^ + 2>2x - 8) , 

^^^(""^ " ^^7 + T^T 

36 [x — 1) 18 [x — 1) 

^ , , X (47a;2 - 79a; + 38) x (-3x3 + 6x - 4) , 

Dh{x) = ^ . - + — 7 — ^logx 

^ ' 108 (a; -1)' 18 (a; -1)^ ^ 

(-43a;2 + lOLx - 52) (2^3 - 9a; + 6) , 

DJx) = ^ o - + - r^logx 

' 108 {x - if 18 {x - if ^ 
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Z°-penguin 



Csm{x) 
Ch{x) 



X (x — 6) X (3x + 2) , 
r H ^ T log X 

X "^1 
8{x-l) ~ %{x-\Y 



16 (y - x) 



logx 7 logy 



8(y-x) 



x-l 

X 



y-1 



x-l 



logx — 



y-1 



logy 



Gluon-penguin 



Esm{x) = 
Eh{x) = 
E^{x) = 



+ llx - 18) (-9x2 + 16a; - 4) , 

+ ^ — —A logx 



12 (x - 1)-^ 6 (x - 1 

X (7x2 _ 292; + 16) a; (3^, _ 2) 



36 (x - 1)'^ 
;-llx2 + 7x-2) 
36 (x - 1)^ 



+ 



+ 



6 (x - 1)' 



logx 



X 



6(x-iy 



logx 
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